The Kondo crossover in shot noise of a single quantum dot 
with orbital degeneracy 
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We investigate out of equilibrium transport through an orbital Kondo system realized in a single 
quantum dot, described by the multiorbital impurity Anderson model. Shot noise and current 
are calculated up to the third order in bias voltage in the particle-hole symmetric case, using the 
renormalized perturbation theory. The derived expressions are asymptotically exact at low energies. 
The resulting Fano factor of the backscattering current Fb is expressed in terms of the Wilson ratio 
R and the orbital degeneracy N as Fb = i+sjjv-ijlfl-i)^ a ^ zero temperature. Then, for small 
Coulomb repulsions U, we calculate the Fano factor exactly up to terms of order U 5 , and also carry 
out the numerical renormalization group calculation for intermediate U in the case of two fold and 
four fold degeneracy (N = 2, 4). As U increases, the charge fluctuation in the dot is suppressed, 
and the Fano factor varies rapidly from the noninteracting value Fb = 1 to the value in the Kondo 
limit Fb = ^fif, near the crossover region U ~ 7rF, with the energy scale of the hybridization T. 

PACS numbers: 71.10.Ay, 71.27.+a, 72.15.Qm 



I. INTRODUCTION 

Kondo physics has been studied since the 1960s for 
dilute magnetic alloys and heavy fermion systems^. Ad- 
vancements in nanofabrication techniques in these years 
have made it possible to experimentally achieve the 
nonequilibrium Kondo state in quantum dots by applying 
a bias voltage, opening a new paradigm of Kondo physics. 
In the early days most research on the nonequilibrium 
Kondo effect has focused on the time averaged current 
or linear conductance^. Recently, the shot noise asso- 
ciated with the Kondo effect has attracted the attention 
of the researchers in this field 4 - - — . 

Shot noise measurements in mesoscopic devices pro- 
vide important information about the effective charge e* 
of current-carrying particles. For instance, the fractional 
charge e* = e/3 of a fractional quantum Hall system 
has been clarified through the shot noise measurement, 
and also the charge e* = 2e of the Cooper-pair has been 
observed in normal metal/superconductor junctions 1 - 1 - - — . 
Furthermore, for other correlated electron systems, the 
shot noise has become an important probe to study the 
properties of the low-energy excitations. In quantum 
dot systems, the observation of a fractional enhancement 
of 5e/3 for the backscattering current in a symmetric 
barrie r 14 ! 15 has stimulated recent studies of the shot noise 
in the Kondo regi me 16 ' 17 . Note that the backscattering 
current is an essential quantity for observing the shot 
noise, and is defined by the deviation from the value of 
the linear current in the unitarity limit. 

Theories for the shot noise in Kondo systems have been 
extended, naturally, to another class of exotic Kondo sys- 
tem : the so-called orbital Kondo effect which has been 
observed in experiments 1 ^—. In the case of multiple 
quantum dots or a single dot with a symmetrical shape, 



the system can have orbital degeneracy. The orbital de- 
grees of freedom have been expected to affect significantly 
the low-energy properties, which can still be described by 
the local Fermi-liquid theory^ 1 -. The shot noise for sys- 
tems with orbital degeneracy has been investigated in the 
Kondo limit, where the local charge degrees of freedom 
in the dot site are quenched^ - — . These studies have 
provided the value of the Fano factor of the backscatter- 
ing current which depends on the orbital degeneracy N, 
and have inspired experimental noise measurements in a 
single-wall carbon nanotube^. 

There have been some qualitative discussions based on 
the Anderson model that capture the physics away from 
the Kondo limit through the local charge degrees of free- 
dom that remain active for finite on-site Coulomb repul- 
sions U— . The approximations used, however, were not 
applicable to low energies, and reliable results in the low- 
temperature Fermi-liquid regime are desired. In a previ- 
ous work, Fujii has derived the expression of the Fano fac- 
tor for the single orbital Anderson model (N = 2) in the 
particle-hole symmetric case using the renormalized per- 
turbation theory (RPT)2&. The RPT is an approach that 
starts with the Fermi-liquid ground state, and gives exact 
asymptotic behavior of the correlation functions at low 
frequencies, low temperatures, and low bias voltages^— . 
Therefore, Fujii's result for the Fano factor i*& is asymp- 
totically exact at low energies, and is expressed in terms 

of a single parameter as Fh = \^\^Ziyi > where R is the 
Wilson ratio that determines the universal Kondo behav- 
ior for all values of the Coulomb repulsion < U < oo^. 

The exact result for the Fano factor for the orbital 
Kondo system with N orbitals has already been given for 
the Kondo limit (U — > oo), as mentioned above. Away 
from the Kondo regime, however, it still has not been 
clarified as to how the Fano factor varies with the value 
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FIG. 1: A schematic plot of the single quantum dot system 
with orbital degrees of freedom. The dot is connected to two 
leads via electron tunneling, conserving the orbital quantum 
number. 



of the Coulomb repulsion. The purpose of the present 
paper is to provide the exact low-energy expression of 
the Fano factor for arbitrary TV and R, on the basis of 
the multiorbital impurity Anderson model. To this end, 
following the derivation in the single orbital case, we use 
the RPT and consider the particle-hole symmetric case^ 8 .. 
Specifically, we start with a general formulation for the 
shot noise, which is based on the nonequilibrium Kubo 
formula for mesoscopic systems^ 3 -, and give the expres- 
sion of the Fano factor in terms of the Wilson ratio R 
and the orbital degeneracy TV. Equations (|28|) and (f30|) 
are the main results of the paper. We also calculate the 
explicit form of the Fano factor for small Coulomb re- 
pulsions up to terms of order C/ 5 , and for intermediate U 
we carry out the numerical renormalization group (NRG) 
calculations for the system with TV = 2 and TV = 4. The 
Fano factor varies rapidly near U ~ ttT from the nonin- 
tcracting value F b = 1 to the value F b = (TV + 8)/ (TV + 4) 
in the Kondo limit, as the crossover from the weak cou- 
pling regime to the Kondo regime takes place. Here, T is 
the linewidth of the dot levels owing to the coupling to 
the leads. 

This paper is organized as follows. In Sec. UH we in- 
troduce a generalized impurity Anderson model for the 
multiorbital quantum dot system, and give a brief expla- 
nation for the RPT to describe the low-energy states. In 
Sec. IIIII wc describe the derivation of the exact expres- 
sion for the current and the shot noise at low energies. 
Then, we discuss the U dependence of the Fano factor 
calculated with the NRG. A brief discussion and sum- 
mary are given in Sec. IIVI 



II. MODEL AND CALCULATION 

A. Multiorbital impurity Anderson model 

Let us consider a single quantum dot system with TV- 
degenerate orbitals. We assume the orbital conservation 
in the tunneling process between the dot and the leads 
in the manner shown in Fig. [1] This assumption has 
been experimentally confirmed to be reasonable for ver- 
tical quantum dot and carbon nanotubc quantum dot 
system a 18 ' 19 . In these assumptions, our system can be 
described by a multiorbital version of the impurity An- 



derson model, 

H = Uo+U!, (1) 

Wo = £ kamcl am c kam + e dm d m d m (2) 

ka.ni m 

+ Yl ( V kccmc{ am d m + h.C.) , 
ham 

Hl = \ E Udld m dl,d m „ (3) 

where d m (dl n ) annihilates (creates) an electron in the dot 
level Cdm with state m, c kam (c\. am ) annihilates (creates) 
a conduction electron with momentum k and state m in 
the lead a = L,R, and U is the on-site Coulomb repulsion 
in the quantum dot. Here, m indicates for the TV-fold or- 
bital state (m = 1, 2, • • • , TV). To investigate the particle- 
hole symmetric Kondo effect, the dot level is chosen to be 
(-dm = — (TV — l)U/2 and the degeneracy TV to be an even 
number. The intrinsic linewidth of the dot levels owing to 
tunnel coupling Vkam is T(w) = (Tl(u>) + Tu(uj))/2 with 
T a (ui) = 2-7T ~ £kam)\Vkam | 2 , which is assumed to 

be independent of the orbital state m. For conduction 
electrons, without any prominent features, T a (oj) does 
not have a strong dependence on w, so it is usual to take 
the case of a wide conduction band with a flat density 
of states limit, where T a (u)) can be taken as a constant 
r Q . It has been shown that an asymmetric barrier gives 
rise to a non-universal shift in the fractional enhance- 
ment for the backscattering current^. To discuss the 
shot noise properties originating from the Kondo correla- 
tion, the lead-dot couplings are assumed to be symmetric: 
r = = r#. The chemical potentials Hl/r — ^ e V/2, 
satisfying fJ*L — (i>r = eV, are measured relative to the 
Fermi level which is defined at zero voltage V = such 
that jjLi, = Hr = 0. 



B. Fermi-liquid description for low-energy states 

We make use of the renormalized perturbation 
theory^, which has been successfully applied to the im- 
purity Anderson model in equilibriu m 3 ^ 40 , and later ex- 
tended to low-bias steady state o 38 ' 42 i 45 . In this section, 
we outline the approach to the multiorbital impurity An- 
derson model described by Eq. (fT]). Furthermore, we pro- 
vide explicit expressions for the renormalized parameters 
for small U obtained exactly up to three-loop contribu- 
tions for general TV. 

The three basic parameters that specify the model are 
the energy level of the dot state e^ m , the linewidth pa- 
rameter r and the Coulomb repulsion U in the dot. The 
low-energy properties of the system can be characterized 
by the quasiparticles with the renormalized parameters, 



3 



defined by 

e d rn = z{t dm + Yr dm {Q)) , (4) 

f = zT , (5) 

U = ^ 2 r^(0,0,0,0) (m^m'), (6) 

where T, dm (uj) is the self-energy of the retarded Green 



function for the dot state: G dm (uj) = [u> — e drn + iT — 
S rfmM] -1 ; z = [1 - 5Sdm( w )/ 9w U=o] _1 is the wave 
function renormalization factor, and T^ m , (oji,u)2, ^3,^i) 
is the local full four-point vertex function for the scatter- 
ing of the electrons with the orbital m and m! . The per- 
turbation theory in powers of U can be reorganized as an 
expansion with respect to the renormalized interaction U, 
taking the free quasiparticle Green's function g r dm {u) = 
(u> — e dm + if) -1 as the zero-order propagator^ - — . The 
three counter terms have to be included to prevent any 
further renormalization of the parameters, e^ m , f, and 
U. Then, the low energy properties can be specified 
entirely in terms of the three renormalized parameters. 
These procedures are the basis of the renormalized per- 
turbation for the impurity Anderson model, presented by 
Hewson^. 

At low bias voltages described by the Fermi-liquid the- 
ory, there are explicit relations between the renormalized 
parameters and the susceptibilities, 
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Sin ( nn drn) 

TV 



f cot(7m rfm ) . 



(7) 

(8) 
(9) 



Here, ridm is the average number of electrons in the dot- 
site, R is the Wilson ratio defined by 



R 



= *d 

' r d 



TV 



(10) 



where x d = 1 - 9S dm( w )/ 9/l L=o 1 / l =o with thc Zeeman 
energy at the dot-site h is the enhancement factor for 
the dot susceptibility, x* c , d = 1 + d Z r drn (u) / de d \ u=0 is 
the enhancement factor for the charge susceptibility, and 
j d = z^ 1 is the enhancement factor for the T-linear spe- 
cific heat coefficient. 

In the particle-hole symmetric case, the electron filling 
is given by n dm = 1/2, and the relations above can be 
simplified as, 



1, 



<~<in 



(11) 



In the Kondo limit (U — > oo), the charge susceptibility 
is suppressed x* c d ~* ano - r can be considered as the 
Kondo temperature as, Tk = ttT/A. Therefore, Eq. (fpQ|) 
shows that the Wilson ratio converges to R — » TV/ (TV — 1) 
in the Kondo limit. In the noninteracting case (U = 0), 
the parameters take the value x* d = X* d — Id — 1 by 
definition. 

Thc behavior of the renormalized parameters for small 
U can be clarified with the perturbation approach. We 
have calculated the enhancement factor "f d and the four- 
point vertex for m m! in the particle-hole symmetric 
case, extending the calculations of Yamada-Yosida for 
TV = 2i&22 to general TV, 
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C(3) - 7 



7T 
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21C(3) 
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^2 



1 - — TV- — +9 
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(12) 



u 4 + 0(u 5 ), (13) 



where u = U/(irT) and £(x) is the Riemann zeta function. Note that 7^ captures the terms of odd order in U for 
TV > 2. Correspondingly, , (0, 0, 0, 0) has finite contributions from thc terms of even order for TV > 2. Thc 
appearance of the zeta function, which is absent in the case of TV = 2, in the coefficients of the perturbation series is 
also caused by thc orbital degeneracy. Substituting Eqs. (|12[) and (fT3j) into Eqs. (0 and (J6j> , and then through Eq. 
(fTTj) . we have obtained the Wilson ratio exactly up to terms of order U 4 , 
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21 



C(3) 
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TV + 77C(3) - — TT 2 -21 



0{u 5 



(14) 
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We note that R is an alternating series up to order U 4 



for TV > 4. 
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For intermediate values of U, the renormalized param- 
eters can be calculated with the NRG approach in the 
case of N — 2 and 4, and the Bcthe ansatz exact solution 
(BAE) in the case of N = 3^^. 



nonequilibrium steady state. The renormalized dot self- 
energy can be calculated exactly up to order cj 2 , V 2 , and 
T 2 , extending the calculations of Ref. |4l| for N = 2 to 
general N. In the particle-hole symmetric case, it takes 
the form, 



III. RENORMALIZED PERTURBATION 
THEORY FOR NONEQUILIBRIUM TRANSPORT 

Here, we derive the expression for the current and 
shot noise of the multiorbital impurity Anderson model 
using the second-order perturbation in the renormal- 
ized interaction U. The obtained expression covers not 
only the strong-couplin g 8 ' 14 ' 15 ' 32 ' 33 and weak-coupling 
limit*^ 7 - but also the whole range of the Coulomb re- 
pulsion U, which is an extention of the previous SU(2) 
resul t 38 ' 41 to general TV. One of the significant advan- 
tages of the RPT is that the current and shot noise cal- 
culated up to the second order in U are asymptotically 
exact at low energies, up to the third order in the bias 
voltage. From the results for the shot noise, we also de- 
rive the exact Fano factor of the backscattering current 
for low bias voltages. 



: (N-1) 
2T 
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0J 2 + -(eV) 2 + {irTf 



(18) 



From this result, T m (uj) can also be determined exactly 
up to order uj 2 , (eV) 2 and T 2 , which leads to the time- 
averaged current at T = up to V 3 from Eq. (fT7|) , 



Ne< 



V 



l + 5(/V-l)(i?-l) 5 
12 



eV 
T 



(19) 



For N — 2, this expression corresponds to Oguri's re- 
sult for the symmetric spin Anderson mode l 41 ' 45 ' 58 . The 
time-averaged current in the Kondo limit also agrees with 
Mora et al.'s result for the half-filled and symmetric lead- 
dot coupling case^ 5 -. 



Current 



The symmetrized current operator J across the quan- 
tum dot system described by Eq. ([3]) is written as, 



J 



Jl - Jr. 



(15) 



where, 

J a. 



- e dt N * 



km 



(16) 



is the current operator for the electrons tunneling from 
lead a to the quantum dot with number operator for elec- 
trons in lead a, N a = £ fcm 4am c w- A general formula 
for the time-averaged current I owing to the applied volt- 
age through a quantum dot has been given by Hershfield 
et al^-£^ and Meir and Wingree n 55 ' 56 . The formula can 
be specialized for the symmetric lead-dot coupling as, 

I= ^ = lJ2 I <^ T ™{") (/iM - /hM) , (17) 



where T m (u>) = — rimG^(cj) is the transmis- 
sion probability for a channel m, and / a (w) = 
[exp((w — fj, a )/T) + 1] 1 is the Fermi distribution func- 
tion for the electrons in lead a. Here, the average (■ • • ) 
is taken over the density matri x 45 ' 57 . In the practical 
calculation, we make use of the density matrix for the 



B. Shot noise 

The current noise S is defined by the correlation func- 
tion for the current fluctuation, 



S= / dt({5J(t),SJ(0)}) , 



(20) 



where the fluctuation operator for a quantity A is given 
by 5 A = A - (A), and {A, B) = AB + BA is anticom- 
mutator. Conventionally, the shot noise has been defined 
by the value of S at zero temperature T = 0, where the 
thermal noise is completely suppressed. This definition of 
the shot noise has been successfully exploited for study- 
ing the properties at low temperatures 5 ^. It is difficult, 
however, for interacting electron systems to separate the 
shot noise and thermal noise at finite temperatures. 

On the basis of a generalized Kubo formalism, Fujii has 
shown that the shot noise can be related to a correlation 
function between the current fluctuation and the charge 
fluctuation, 



.5',, 



({SJ,e(SN L -SN R )}) , 



(21) 



and has suggested this correlation function as a finite- 
temperature shot noise.—. The correlation function Sh 
has been shown to satisfy the identity, 



S h = S-Ak B TG 



(22) 



with the differential conductance G and S defined in Eq. 
([20| . The expression (|22|) is identical to an empirical 
formula that has been applied as an estimation of the 
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effective shot noise at finite temperatures from the mea- 
surement of the current noise and the differential conduc- 
tance. 

It has also been confirmed that a number of the proper- 
ties of the shot noise can be rederived from the expression 
of Sh in Eq. ([21]) . First, at zero temperature, Eq. ([22]) is 
simplified as Sh = S. Therefore, Sh defined in Eq. (f2"Tj) 
clearly agrees with the conventional shot noise defined at 
T = by S. Second, it has been proved that the shot 
noise defined by Eq. (|2ip vanishes in the linear response 
regime (V = 0): Sh = . This means that the Nyquist- 
Johnson relation is also reproduced. These observations 
show that Eq. ([21]) can be regarded as an extension of 
the shot noise to all temperatures. 

Therefore Sh given in Eq. (|2"T]) enables one to study 
finite-temperature effects on the shot noise directly, with- 
out calculating S and G separately In the present pa- 
per, however, we focus on the shot noise at zero tem- 
perature. Specifically, we calculate Sh exactly up to V 3 
in the particle-hole symmetric case where the occupa- 
tion number of each orbital is given by n dm — 1/2. We 
demonstrate two different strategies to achieve our goal. 



The first starts from Sh defined in Eq. ([21]) and reaches 
the result given in Eq. ([2"5]) . The second from Eq. ([2"0]) 
yields the final expression given in Eq. ([2"5)> . The results 
obtained in these two ways, Eq. ([2"5]) and Eq. ([2"5f. agree 
with each other naturally, although the contributions of 
each Feynman diagram for Sh and that for S do not have 
one-to-one correspondence as summarized in Table [IJ 

1. Shot noise Sh at absolute zero 

We now calculate the shot noise from the expression 
given in Eq. (|21[) and derive the leading asymptotic de- 
pendence of applied bias voltage at T = in the particle- 
hole symmetric case. 

First, substituting the current operator Eq. (|15[) and 
charge fluctuation operator into Eq. ([2"T]). the shot noise 
is readily expanded in the Keldysh formalism as, 

S h = [(FhLL - F hLR ) + (L^ R)] + (c.c.) , (23) 

with, 



= -*T E ^ m (r c 5 c cL rn (0+)d m (0 + )4., Q , m ,(0-) Cfe , Q w(0-)) (24) 

li — \ / connected 

k,k' ,m,m' 



Here, S c = T c exp[— i J c dt Hi(t)] is the time-evolution op- 
erator, the Keldysh contour c runs along the forward time 
direction on the branch "— " followed by the backward 
evolution on the branch "+" , T c is the corresponding 
contour-ordering operator, and (. . . ) conncctoc i takes the 
sum of all connected diagrams. The equal time correla- 
tion is defined by the Keldysh contour: in Eq. ([2"4")l is 
on the branch '±'. The form of Eq. ([M|) enables us to ap- 
ply perturbation expansion in U. Then, the asymptotic 
behavior of the shot noise up to V 3 can be calculated 
using the RPT. 

The contributions of the second order perturbation in 
U can be classified using the diagrams shown in Fig. [5] 
Note that for Sh defined in Eq. (|21[) the fluctuation op- 
erator assigned for the left end of each diagram and that 
for the right end are different, i.e., one of the two is the 
current operator and the other is the charge fluctuation 
operator. The calculated contributions of each diagram 
are summarized in Table[U In the particle-hole symmetry 
case the contributions from diagram (d) and (e) vanish. 
The contribution from diagram (g) and that from (f) are 
the same for arbitrary parameters. Furthermore the con- 
tribution from (h) coincides with that of (f ) and (g) in the 
particle-hole symmetric case for the symmetric lead-dot 
coupling. 




(f) (9) (h) 




FIG. 2: The Feynman diagrams for Sh and S generated from 
Eqs. (|24|l and (|27|) . The contributions up to second order in 
U are shown. The Hartree term is included in the nonpertur- 
bative Green's function. Note that the fluctuation operators 
of the two verticies at the left and right ends for Sh and those 
for S are different. 

shot noise at T = up to order V 3 , 

s h = + 4 b) + 4 c) + 3^ 
-^ii(£)"[^-w-4 



Collecting all these contributions, we obtain the total 



(25) 
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TABLE I: The contributions of each diagram shown in Fig.[2]for S h and S. The unit is C = (2e 3 /h)\V\(eV/T) 2 . 



diagram 


(a) 




(b) 






(c) 


(d),(e) 




(f),(g),(h) 


s h 


NC/12 


N(N — 


1)(R- 


l)' z C/2A 


5N{N - 


i)(r- iyc/24 





N(N 


- l)(R-lfC/Q 


s 


NC/12 


N(N - 


-1)(R- 


- l) 2 C/8 


N(N — 


1)(R- l) 2 C/8 





N(N 


- 1){R- l) 2 C/6 



For N = 2, this result agrees with the one for the SU(2) 
Anderson mode l 38 ' 60 . In the Kondo limit where R — > 
N/(N — 1), our result agrees with Mora et al.'s result in 
the particle-hole symmetric case for the symmetric lead- 
dot coupling 3 ^. 

2. Current noise S at absolute zero 

Here we calculate the shot noise via the current noise 
S denned in Eq. (|20|) . Because the thermal noise is com- 
pletely suppressed at zero temperature, we can extract 



the pure shot noise from the current noise S. We derive 
the asymptotic form of the shot noise along a similar line, 
as that described above. 

Substituting the current operator Eq. (|15p into Eq. 
(f2"0)) readily leads to the current noise expression in the 
Keldysh formalism as, 

+ (c.c.) , (26) 

with, 



L \ / connected 

-V ak V*, k , (T c S c ci km (t)d m (t)di,(t>)c a , k , m ,(t>) 
I 



connected 



(27) 



The contributions of the second order perturbation in U 
can be classified using the diagrams shown in Fig. [2] for 
this correlation function. Note that for the current noise 
S both of the two verticies, at the left and right ends of 
each diagram, are given by the current fluctuation oper- 
ator 5 J. The calculated contributions of each diagram 
are summarized in Table HI For this reason, the contri- 
butions from the diagram (b) and those from (c) are the 
same for S, in contrast to the case for Sh- Nevertheless, 
there are some similarities between the two cases. The 
contributions from (d) and (e) vanish in the particle-hole 
symmetric case. The contribution from diagram (g) and 
that from (f) are the same for arbitrary parameters, and 
the contribution from (h) coincides with those of (g) or 
(f) in the particle- hole symmetric case for the symmetric 
lead-dot coupling. In the Kondo limit the contribution 
of each diagram agrees with that of the corresponding 
diagram used in Mora's calculations 3 ^. 

Collecting all these contributions, we obtain the total 
current noise up to order V 3 as, 



S = 5 (a) + 2S {b) + 3S (f) 



27Ve 3 /eV 

h 1 It 



^ + ^(7V-l)(i?-l) 2 



(28) 



This result agrees with Eq. (|25p which has been deduced 



from Sh- 



Fano factor 



We consider the backscattering current lb which con- 
tains all effects of the quantum and thermal fluctuations. 
It is denned via the deviation of the nonequilibrium cur- 
rent / from the value in the unitary limit Ne 2 V/h, as, 



h 



Ne 1 
Ne 2 



V -I 

-v eV 



1 

12 



5(iV- 



12 



(29) 



Here, we have used the result for / given in Eq. (|19[) . 
Generally, the current noise for the forward current and 
that for the backscattering current are equivalent in the 
case where the systems is coupled to two leads. There- 
fore, the Fano factor of the backscattering current can be 
expressed in the form, 



F h = 



S i + 9(jy- l)(fl-l) 2 
2el b ~ l + 5(A^-l)(i?-l) 2 



(30) 



and this is one of the main results of the present work. 

In the nonintcracting case (U = 0), the Wilson ratio 
takes the value R = 1, which leads to the Poisson noise 
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of S = 2elb and Ff, = 1. For small {/, substituting the 
perturbation expansion for the Wilson ratio given in Eq. 



J 



(jl4| into Eq. (|30|) . the Fano factor can be calculated ex- 
actly up to terms of order U 5 , 



F b = l + 4(/V -l)u 2 -8(N -l)(N -2)u 3 + 2{N - 1) 
-(N- l)(iV-2) 



67V 2 - (26 + tt 2 )N - 3tt 2 + 66 



178 



167V 2 + (84C(3) - IOtt 2 - 160) N + 616C(3) - — tt 2 + 8 



it 5 + 0(w 6 ) 



(31) 



r 



TABLE II: The Fano factor F b in the Kondo limit (U 
for several choices of orbital degeneracy N. 



N 


2 


4 


6 


8 


— > oo 


F b 


5/3 


3/2 


7/5 


4/3 


-)■ 1 



This expression obviously shows that a larger orbital de- 
generacy makes the initial rise of Ft steeper, as the coef- 
ficient of the order u 2 term increases with N. 

In the opposite limit, namely, the Kondo limit (U — > 
oo), the Wilson ratio approaches to the universal value 
R — > N/(N — 1), and the Fano factor takes the form, 



F, 



N + 8 
7V + 4 



(32) 



The explicit values for several N are given in the Table HU 
Specifically in the limit of large orbital degeneracy N — > 
oo, the shot noises take the Poisson value S — s- 2elb, and 
the Fano factor approaches to the nonintcracting value 
F(, — > 1 even though the Coulomb repulsion has been 
taken first to be U — >• oo. This originates from the fact 
that renormalization is weakened in the limit of large TV. 
In practice, renormalized parameters takes ?7 — >■ and 
F -> T in the large N limit (R -> 1) of Eqs. and 
©. For N = 2, i.e., in the SU(2) case, our result is 
consistent with the previous results obtained by Gogolin 
and Komnik, Sela et al. and Fujii.— 

In order to clarify the behavior of Fb in the intermedi- 
ate values of U, we have carried out the NRG calculations 
for N — 2 and 4. Specifically, we have deduced the Wil- 
son ratio R from the low-energy NRG fixed points. This 
method has also been applied to the two-channel Ander- 
son model in the recent work of Nishikawa el al^. In 
the case of N = 2, the NRG and Bethe ansatz results for 
the renormalized parameters have been shown to agree 
wel l 52 i 53 ' 63 as seen in Fig. [3] (a). There is no Bethe ansatz 
solutions for N > 2 in the particle-hole symmetric case, 
still the NRG is applicable for N = 4. We can see, in 
Fig. [3] (a), that the NRG results for the Wilson ratio for 
N = 4 are in good agreement with the perturbation re- 
sult up to a term of order C7 4 given in Eq. (fl"4| for a small 
Coulomb repulsion u < 0.3. Furthermore, the NRG re- 
sults approach the correct universal value R — > 4/3 for 
larger U . Therefore, the U dependence of the Fano factor 



for N = 4 can also be deduced, using expression (|3U)) with 
the nonperturbative NRG approach, which is discussed 
below. 

As the Coulomb repulsion increases further, the charge 
fluctuation in the dot is suppressed and the Wilson ratio 
increases rapidly from the nonintcracting value R = 1 to 
the universal value in the Kondo limit. This crossover 
from the weak-coupling regime to the Kondo regime is 
also observed in the Fano factor, shown in Fig. [3] (b), as 
a rapid convergence of the Fano factor to the universal 
value Fb = 5/3 and 3/2 for N = 2 and 4, respectively. 
Note that the Wilson ratio is a decreasing function of 
X* Jx*di namely, the ratio of the charge susceptibility to 
the susceptibility as shown in Eq. (fTU]) . Therefore, the 
charge fluctuation at the dot-site suppresses the Fano 
factor for small Coulomb repulsion U < 71T. The NRG 
results for N — 4 have successfully revealed the precise 
feature of the crossover for the system with the orbital 
degeneracy 

The perturbation series is applicable also for large 
orbital degeneracies N where the NRG is not feasible. 
Therefore, it is worthwhile to discuss the convergence of 
the series expansion. For this purpose, the Fano factor 
for N = 4 is plotted in Fig. 0] (a), keeping the first few 
terms of the series given in Eq. (f3"T|) . We can see that the 
values of Fb deduced from Eq. (|3"Tj) agree with the NRG 
result for a relatively narrow range of the Coulomb inter- 
action u < 0.1, while the series expansion for the Wilson 
ratio quantitatively works in a wider range u < 0.3. This 
discrepancy is caused by the expansion of the denomi- 
nator on the right-hand side of Eq. (|3"U| with respect to 
(iJ — 1), the convergence radius for which becomes small 
as shown in the Appendix. This can be resolved, how- 
ever, using the expression of the Fano factor Eq. (f30]) . 
as it is, without expanding the denominator, and then 
substituting there the series expansion for the Wilson ra- 
tion given in Eq. (|14[) . The results obtained in this way 
agree with the NRG data for u < 0.3, as shown by the 
dashed line in Fig. @] (a). Specifically, for N = 4, both 
the coefficient of the order u 3 term and that for the or- 
der w 4 term in Eq. (f3"Tj) become negative, while that for 
the order u 5 term is positive. For this reason, in Fig. [4] 
(a), the fourth order curve deviates from the NRG re- 
sult earlier than the third order curve, and then the fifth 
order curve again approaches the NRG result. Finally, 
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FIG. 3: (Color online) (a) The Wilson ratio R, and (b) the 
Fano factor F\> as functions of the on-site Coulomb repulsion 
U . The dashed-dotted line denotes the results deduced from 
the BAE for N = 2. The NRG results are plotted for TV = 2 
(a), and for TV — 4 with the solid circle (•). The dashed 
line in (a) is the Wilson ratio for TV = 4, obtained from the 
perturbative expansion up to terms of order U 4 given in Eq. 



we examine the TV = 6 case, for which nonperturbative 
approaches are not available at present. The value of the 
Fano factor for TV = 6 deduced from the perturbation 
series are shown in Fig. [4] (b) . We see that the results 
are quantitatively valid for u < 0.1, although the conver- 
gence of the series expansion for R and that for Fb given 
in Eqs. (fT4"|) and (piTj) become worse with an increase of 
TV. Nevertheless, the results show clearly that the initial 
rise of Ft near u ~ 0.1 is steeper for TV = 6 than that for 
TV = 4. 



IV. DISCUSSION AND SUMMARY 

At the end of the previous section, we have calculated 
the value of the Fano factor from the data for the Wilson 
ratio. Conversely, it may be possible to determine the 
value of the Wilson ratio from the measurements of the 
Fano factor, using the expression, 

« = 1 + l/ ("-w»W (33) 
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0.05 0.1 0.15 0.2 



U/ttT 

FIG. 4: (Color online) Perturbation results for the Fano factor 
for (a) TV = 4, and (b) TV = 6. The dashed line is obtained just 
by substituting Eq. (TT4l into Eq. (TSUI) . The dotted, dashed- 
dotted, and dashed-dot-dotted lines denote the values from 
the series expansion up to terms of order U 3 , U 4 and U 5 , 
respectively, given in Eq. (|31[) . The solid line with the solid 
circle (•) is the NRG results for TV = 4. 



which is readily obtained by inverting Eq. (f30|) . This 
scheme for estimating R from Ff, does not need the value 
of the g-factor for the electrons in the dot. Only the ob- 
servation of Fb determines the accuracy of observed R. 
Alternatively, one can estimate the Wilson ratio from the 
measurements of the susceptibilities x* d/Xd- ^ n ^ n ^ s case j 
the g factor is also needed to determine them accurately. 
Therefore, Eq. may enable us to experimentally ob- 
serve the Wilson ratio with a more reasonable accuracy. 

Finally, we comment on a physical interpretation for 
the values of the Fano factor for the orbital Kondo sys- 
tems. The expression for the Fano factor presented in 
Eq. (|30[) shows that effects of the Coulomb repulsion 
enter through the Wilson ratio R, which characterizes 
the low-energy Fermi liquid state. The terms having R 
in the coefficient come from the vertex corrections and 
self-energy appearing in the RPT approach to the shot 
noise defined by Eqs. (f20|) or (|2ip . and the nonequilib- 
rium current /. In this aspect, the role of the Fano factor 
for the shot noise resembles the role of the Stoner factor 
for the magnetic susceptibility. These factors quantita- 
tively enhance the response of the system against exter- 
nal fields, but do not change the qualitative feature of 
the low-energy states. Therefore, the deviations of the 
Fano factor from the noninteracting value are caused by 
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the scattering between the renormalizcd quasiparticles of 
the Fermi liquid, rather than an effective charge as with 
other kinds of elementary excitations. 

In summary, we have studied the current and shot 
noise through the single quantum dot described by 
the multiorbital impurity Anderson model, using the 
nonequilibrium Kubo formalism. Employing the renor- 
malizcd perturbation theory, the current and the shot 
noise have been calculated exactly in the particle-hole 
symmetric case up to order V 3 at T = 0. The result for 
the shot noise is given by Eq. (|2"5)) , and then the Fano fac- 
tor Ft of the backscattering current is determined by two 
parameters, i.e., the Wilson ratio R and the orbital de- 
generacy N, as shown in Eq. (|30|) . We have also presented 
the explicit form of the Fano factor for small Coulomb 
repulsions up to terms of order U 5 . Furthermore we have 
deduced the dependence of Ft for N = 2 and A" = 4 on 
the Coulomb repulsion, using the Wilson ratio obtained 
with the NRG approach. The Fano factor varies mono- 
tonically with an increase of U from the noninteracting 
value Ft = 1 until almost saturating for U > irT to the 
value in the Kondo limit, F b = (N + 8)/(A + 4). For 
small U, the charge fluctuation at the dot-site suppresses 
the Fano factor. 
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series in U given in Eq. (fT4|) slowly converges for large 
N. Therefore, it can be concluded that the convergence 
of the Fano factor as a power series in U given in Eq. 
(f3"T|) is slower for larger degeneracy N. 

Second, we consider the radius of convergence of Eq. 
([31]) . As mentioned in the text, the value of the Wilson 
ratio is bounded in a range, 

0<R-l<-±_ 7 (A3) 

for < U < oo. For A" = 2,4, it follows from Eqs. 
(|A2[) and (|A3|) that the radius of convergence of the Fano 
factor as a power series of U is determined by, 

^x-^Jsuhr)- (A4) 

Therefore, the radius of convergence can be obtained as 
u ~ 0.48 for N = 2 and u ~ 0.72 for N = 4 from 
the Wilson ratio data shown in Fig. [3] (a). However, 
for A" > 6, the range of R — 1 given in Eq. (|A3j) is 
narrower than that of Eq. (|A2[) . Therefore, the Fano 
factor for A^ > 6 can be written as a power series of 
U for arbitrarily positive U, similarly to the universal 
quantities of the Kondo effect. 
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Appendix A: The convergence of the Fano factor 

m 

We discuss the convergence of the Fano factor as a 
power series of U given in (|31|) . 

First, the Fano factor given in (|5D| can be expanded 
as a power series in R — 1 around R — 1 = 0, as, 

oo 

F b =[l + 9(N - 1)(R - l) 2 ] Cn(N)(R - l) 2n ,(A1) 

71=0 

with C n (N) = [—5 (N — 1)]™. and the radius of conver- 
gence, 

v-^Jjvh)- (A2) 

It follows from Eq. (|A1[) that for large N the power series 
slowly converges. In addition, the Wilson ratio as a power 
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